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, The exclusive processes '^H{e, e'p)n, ^He{e, e'p)'^H and ^He{e, e'p){pn), have been analyzed using 

^ ' realistic few-body wave functions and treating the final state interaction (FSI) within a Generalized 

Eikonal Approximation (GEA), based upon the direct calculation of the Feynman diagrams de- 
scribing the rescattering of the struck nucleon with the nucleons of the A — 1 system. The approach 
represents an improvement of the conventional Glauber approach (GA), in that it allows one to 



> 
O 

' take into account the effects of the nuclear excitation of the A — 1 system on the rescattering of the 

■ struck nucleon. Using realistic three-body wave functions corresponding to the AFIS interaction, 
the results of our parameter free calculations are compared with available experimental data. It is 

Q I found that in some kinematical conditions FSI effects represent small corrections, whereas in other 

Ch ' kinematics conditions they are very large and absolutely necessary to provide a satisfactory agree- 

, ment between theoretical calculations and experimental data. It is shown that in the kinematics 

X, 

^ , of the experimental data which have been considered, covering the region of missing momentum 

and energy with p^ < 0.6 GeV/c and < 100 MeV in the perpendicular kinematics, the GA 
and GEA predictions differ only by less than ~ 3 — 4%. 
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I. INTRODUCTION 



One of the main aims of nowadays hadronic physics is the investigation of the hmits of 
vahdity of the so called Standard Model of nuclei, i.e. the description of nuclei in terms of 
the solution of the non relativistic Schrodinger equation containing realistic nucleon-nucleon 
interactions. To this end, exclusive lepton scattering could be very useful, for it might 
yield relevant information on the nuclear wave function, provided the initial and final states 
involved in the scattering process are described within a consistent, reliable approach. In 
the case of few-body systems, a consistent treatment of initial and final states is nowadays 
possible at low energies (see e.g. and References therein quoted), but at higher energies, 
when the number of partial waves sharply increases and nucleon-nucleon (NN) interaction 
becomes highly inelastic, the Schrodinger approach becomes impractical and other methods 
have to be employed. In the case of complex nuclei, additional difficulties arise due to the 
approximations which are still necessary to solve the many-body problem. As a matter 
of fact, whereas fundamental progress has been made in recent years in the calculation of 
various properties of light nuclei (see e.g. 0| and Q] and References therein quoted), much 
remains to be done for the treatment of the continuum, for which various approximate 
treatments of the final state cannot be avoided. In this context, it should be stressed that 
calculations involving few-body systems, where the ground state can be treated exactly, can 
also be very useful to investigate the limits of validity of various approximate schemes to 
treat the continuum and their possible extension to complex nuclei. 

The aim of this paper is to present the results of a systematic theoretical investigation 
of the exclusive process A{e, e'p)B off (to be also denoted by D) and ^He, based on a 
reliable description of: 

1. initial state correlations (ISC), treated by the use of the status-of-the-art few-body 
wave functions 0| corresponding to the yiyi8 interaction Qj; 

2. final state interactions (FSI), treated within a relativistic framework based upon the 
calculation of the relevant Feynman diagrams which describe the rescattering of the 
struck nucleon by the other A — 1 spectator nucleons of the target. 

Whereas a correct treatment of ISC in few-body systems is automatically achieved by 
the use of realistic wave functions, the treatment of FSI at high energies is still matter of 
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discussions. The approach we are going to use has several non trivial advantages, in that it 
allows one to work within a relativistic framework provided by the use of Feynman diagrams 
and, moreover, it can be applied, in principle, to the treatment of exclusive A{e,e'p)B pro- 
cesses off complex nuclei as well. It should be stressed, at this point, that the diagrammatic 



approach we are talking about is not a new one: it has been first formulated in Ref. pi and 



(see also Ref. 



within a spin-less treatment of particle-nucleus scattering, and applied 



subsequently to various types of high energy processes with nuclear targets. More recently, 
the diagrammatic approach has been generalized to the treatment of the FSI in exclusive 
, Inl and Me, e'2p]B 



A(e, e'p)B 



also been used in Ref. 



he; 



processes, and a Feynman diagram approach has 



and 1^ to take into account off-shell effects both in inclusive, 
A{e, e')X, and exclusive, A(e, e'p)B, processes. 

The diagrammatic ap pro ach we are referring to, is a generalization of the standard 
Glauber Approach (GA) [l6| based on the eikonal approximation, so that, following Refs. 
[l^ . we will call it Generalized Eikonal Approximation (GEA). 

It is well known that the application of the GA to the treatment of A{e, e'p)B processes 
requires the following approximations: i) the NN scattering amplitude is obtained within 
the eikonal approximation; ii) the nucleons of the spectator system A — 1 are stationary 
during the multiple scattering with the struck nucleon (the frozen approximation) , and 
iii) only perpendicular momentum transfer components in the NN scattering amplitude are 
considered. In the GEA the frozen approximation is partly removed by taking into account 
the excitation energy of the A — 1 system, which results in a correction term to the standard 
profile function of GA, leading to an additional contribution to the longitudinal component 
of the missing momentum. 

In the present paper we apply both the GA and the GEA to the calculation of the 
processes '^H{e,e'p)n, ^He(e,e'pYH, and ^Heje, e'p)(np), and compare our results with 



Q Q, Q 0, 



The ^He wave function of the Pisa 



available experimental data 
group [2] , corresponding to the AV18 interaction , will be used in the calculations. We will 
not consider, for the time being. Meson Exchange Currents (MEG), A-Isobar Configurations, 
and similar effects, which have been the object of ^intensive theoretical studies in A (e, e 'p)B 
processes off both few-body systems (see e^. 
References therein quoted). As in Refs 



. 23^2' 

m 



24 n and complex nuclei (see e.g. j25[ and 



12j |. we fully concentrate on the effects 



of the FSI but we will consider kinematical conditions for which the effects from meson 
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exchange currents (MEC) and A excitation effects are expected to be small corrections, 
and, whenever possible, we will compare our results with the results by other authors which 
include these effects. 

The structure of the paper is as follows: in Section El the basic formalism of lepton- 
hadron scattering is briefly illustrated and the main formulae are obtained; in Section UTTl the 
concept of Plane Wave Impulse Approximation (PWIA) and Spectral Function are recalled; 
in Section IV, the GEA is introduced, the relevant Feynman diagrams which one needs to 
take into account in the treatment of the full FSI are analyzed, and the problem of the 
factorization of the lepton-nucleus cross section within the GA and GEA is also discussed; 
the results of the calculations and their comparison with available experimental data are 
shown in Section |3 eventually, the Summary and Conclusions are presented in Section IVTl 
Some details concerning the formal aspects of our approach are given in Appendices A and B. 
Preliminary results of our calculations have been reported in Ref. j3| and |27|. 



II. BASIC FORMULAE OF (e, e'p) SCATTERING OFF NUCLEI 

The one-photon-exchange diagram for the process A{e, e'p){A — 1), where A — 1 denotes 
a system of A — 1 nucleons in a bound or continuum state, is presented in Fig. Q where the 
relevant four-momenta in the scattering processes are shown, namely the electron momenta 
before and after interaction, k = (-E, k) and k' = (i?',k'), the momentum of the target 
nucleus Pa = {Ea, Pa) and, eventually, the momenta of the final proton and the final A — 1 



system, pi = {^pi^ + M^,pi) and Pa-i = (^Pa-i + Pa-i), where Mn is the 

nucleon mass, Mj4_^ = Ma-i + E^_^, and -E^-i the intrinsic excitation energy of the 
A—1 system. 

Let us briefly recall some useful formulae regarding the process described by the diagram 



shown in Fig. 
form (see e.g. 



The differential cross section for the exclusive process has the following 

1) 

(fa 

E'^ cos2 I 

where a... ^ .s the Mott cross section, a the fine-structure constant, Q ^ 

—q^ = —{k — k'Y = — Q'o = '^EE' sin^ 9 /2 the four-momentum transfer, 9 = 9^, the 
scattering angle. The quantities Ifj^i, and are the reduced leptonic and hadronic tensors. 



respectively; the former has the well known standard form ( 28|), whereas the latter can be 
written as follows 

AttMa „^ 

X {oaP A\J^{0)\aNPi, aA-iP a-iE{_i) {E{_-^P A-iaA-i,PiaN\Ju{0)\(yAP a) , (2) 

where denotes the set of discrete quantum numbers of systems A, A — 1 and A^. In 
Eq. (j21) the vector {a^pi, aA-iP a-iEa^i) consists asymptotically of a nucleus A — 1, with 
momentum Pa-i and intrinsic excitation energy Ea-i, and a nucleon with momentum pi. 
Two relevant experimentally measurable quantities which characterize the process are the 
missing momentum (i.e. the momentum of the A — 1 system), and the missing energy 
Em defined, respectively, by 



Pr 



q - Pi Em = JPLi + Mn-Ma = Mn + Ma-i - Ma + E{_, = E^in + ^Li, (3) 



where Emin = Ea — Ea-i = Mjq + Ma-i — Ma, and the (positive) ground-state energies 
of A and A — 1 are denoted by Ea and Ea~i, respectively. The exclusive cross section can 
then be written in the well-known form 

(fa 

where i = {L,T, LT,TT}, and Vl, Vt, Vlt, and Vtt are well-known kinematical factors. 

The evaluation of the nuclear response functions W/^ requires the knowledge of the nu- 
clear vectors I^aPa) and \aNPi, o:a~iP a~iEa-i) and the nuclear current operators J^(0). 
Nowadays, there is no rigorous quantum field theory to describe, from first principles, a many 
body hadronic system, and one is forced to adhere to various approximations. Whereas at 
relatively low energies a consistent non relativistic treatment of the electro-disintegration of 
two- and three-body systems can be pursued, with increasing energy the treatment of the 
three- body final state requires proper approximations. In the present paper we describe the 
two- and three-body ground states in terms of realistic wave functions generated by modern 
two-body interactions [^, and treat the final state interaction by a diagrammatic approach 
of the elastic rescattering of the struck nucleon with the nucleons of the A — 1 system. The 
relevant diagrams which, within such an approximation, replace the One-Photon-Exchange 
diagram of Fig. 1, are shown in Fig. |21 the first one represents the Plane Wave Impulse 
Approximation (PWIA), whereas the other ones the final state rescattering (FSI). 



Although the PWIA appears to have a hmited range of vahdity, it is useful to analyze 
its predictions since, within such an approximation, the cross section is directly related to 
a quantity, the Spectral Function, which, in the case of few-body systems, can be calculated 



with high degrees of accuracy (see |29l. l30L l3ll.l32j). The relevant point here is that, provided 



the FSI of the struck nucleon with the A—1 system can be disregarded, the Spectral Function 
yields direct information on the nuclear wave function. For such a reason, we will present 
our results obtained within two distinct approaches: 

1. the PWIA (Fig.|2K)), when the struck proton is described by a plane wave, whereas the 
A — 1 system in the final state, with momentum Pa-i, represents the bound or continuum 
state solutions of the Schrodinger equation with the same potential used to obtain the A- 
body wave function (note that some authors call PWIA the state in which all particles in 
the continuum are described by plane waves); 

2. the full FSI approach (Fig. |2b,c)), when the A — 1 system (in the ground or continuum 
states) is still described by the exact solution of the Schrodinger equation, and the interaction 
of the struck nucleon with the A — 1 nucleons is treated by evaluating the Feynman diagrams 
of Fig. 2, either in the GA or the GEA approximations. 



III. THE PLANE WAVE IMPULSE APPROXIMATION AND THE NUCLEAR 
SPECTRAL FUNCTION 

The main merit of the PWIA is that it allows one to express the nuclear response functions 
W/^ in terms of the nucleon response functions which are very well known from e — N 
experiments. As a matter of fact, by expressing the hadronic tensor for the nucleus A (Eq. 
(121)) in terms of the hadronic tensor for the nucleon 

= 7-l:rEE(2vr)V^)(p + g-Pi)(«;vPi|J,^(o)|a;vP;)(«;vP;i^.^ (5) 



N 



the cross section assumes the following form (see e.g. Refs |33L l34l l35l|) 

;^^^ = if(Q^..p„).-(Q^p,„)p.(|k,l.s). (6) 

where =q' - ql ( = * + Ma - ^J{H + (M{_iY - y'k,^ + ), and K(Q\ x. p) a 
kinematical factor. In Eq. a^^ {Q'^ , pm) is the cross section describing electron scattering 
by an off-shell nucleon, x = Q^/2MNqQ is the Bjorken scaling variable, ki = —pm is the 



nucleoli momentum before interaction, E = Emis = Emin + -E'a-i removal energy and 

P(|ki|,ii^) is the nucleon Spectral Function, which can be written as follows: 



{2Tif2JA + 
X 6(E- {E{_, + E, 



Y E E 

/ Ma,Ma~i,'^n 



(aAPA|ttjvki,aA--iPA-i^i_i) X (7) 



where Ma, M.a-i, and ct„, are the spin projections, and the sum over / includes all possible 
discrete and continuum states of the A ~ 1 system. 

Whereas the Spectral Function for the Deuteron (D) has a particularly simple form, viz. 



Pi,(|ki|,E) =nz)(|ki|)5(E-eB), 



(8) 



where en is the (positive) binding energy of the deuteron and nD(|ki|) = 
(27r^) ^ (u|(|ki|) +M£5(|ki|)) the nucleon momentum distribution, in the case of A=3, the 
proton Spectral Function consists of two parts. 



PHe(|kl|,S) = P,.(|ki|,E)+P,,.(|ki|,E), 

The first one, or ground part Pgr, has the following form 

Pgr{\kllE) = ngr{\kl\)5{E - E^in ) , 



(9) 



(10) 



where Emin = {E^l — \E2\ ~ 5A9MeV, and n^rdkil), which corresponds to the two-body 
break-up (2bbu) channel ^He ^ D + p, is (hereafter, the projection of the spin of nucleon i 
will be denoted by Si) 

^..(|ki|) = 7^^ E [e-^P'^^xl^^^'\r)^^:{p,r)dpdr\ (11) 

In Eq. (HH) ^^/(P, r) is the ^He wave function, M.^ the projection of the spin of '^He , and 
r and p the Jacobi coordinates describing, respectively, the motion of the spectator pair and 
the motion of the struck (active) nucleon with respect to the CM of the pair. 

The second, or excited part Pe^,., of P//e(|ki|, E), corresponds to the three-body break-up 
(3bbu) channel ^He —>■ {up) + p and can be written as follows 

I I ^ r dH 



Pex(|kl|,P) 



X 5\E 



M, 



(27r)3 



e^'P'''x\,yi{v)K!{P,^)dpdv 



N 



(12) 



where ^*p(r) is the two-body continuum wave function characterized by spin projection 5*23 
and by the relative momentum t = of the np pair in the continuum. Obviously, for 

the neutron Spectral Function, only the excited part (fT^ contributes. 

In Fig.iniwe show the Spectral Function of '^He obtained using the variational three-body 
wave function by the Pisa group 2] corresponding to the realistic AV18 potential \^ (see 
Appendix A). The two-body wave function entering Eq. ()12j) has been obtained by solving 
the Schrodinger equation for the continuum using the same Al^lS two nucleon potential. Our 
results for the Spectral Function agree with the ones obtained in Ref. jsil, where the same 
three-body wave function has been used. The normalization of the Proton Spectral Function 
has been fixed to 2 (two protons) and the normalization of the neutron Spectral Function 
to one. In Fig. IHl we also show the results predicted by the Plane Wave Approximation 
(PWA), which corresponds to the replacement of the continuum interacting {n — p) pair 
wave function with two plane waves. The left panel of Fig. El shows Ug^ and the right panel 
Pex- It can be seen that: i) Pex exhibits maxima centered approximately at Em 



ii) around these values of Em and ki the Spectral Functions, calculated disregarding the 
interaction in the A^A^-pair in the continuum (PWA) and taking it into account (PWIA), 
are almost identical, in agreement with the results obtained l ong ago 3^ with the Spectral 
Function corresponding to the Reid Soft Core Interaction 3y|. The region centered at 
Em ~ \. is the so-called two-nucleon correlation region when one of the nucleons 

of the spectator A^A^-pair is fast, the other one being basically at rest ( for an improved 
description which takes into account the motion of the third, uncorrelated nucleon, or the 
A — 2 spectator system in case of heavier nuclei, see 38|). Then the fast nucleon becomes 
strongly correlated with the active nucleon (the proton, in the case of the proton Spectral 
Function, or the neutron, in the case of the neutron Spectral Function) forming a correlated 
pair which carries most of the nuclear momentum. In this case, it is intuitively expected 
that the slow nucleon acts as a passive spectator and, consequently, only the interaction 
in the correlated pair can be relevant for the Spectral Function. Hence, in this region the 
calculations including or omitting the interaction in the spectator pair, are expected to 
provide essentially the same results, as confirmed by present and previous calculations of 
the Spectral Function j3, 22,1^. The situation which has been just described, is clearly 
illustrated in Fig. U where the three-dimensional neutron Spectral Function is presented. 
The PWIA results suggest that experimental insight about the structure of the nuclear 
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wave function at short distances can be obtained from A{e,e'p){A — 1) processes, provided 
the PWIA entirely exhausts the reaction mechanism. Unfortunately, we know that in many 
cases the simplified PWIA mechanism fails to describe the experimental data (see, e.g. a 
recent discussion in Ref. However, a properly chosen kinematics can still leave room 

for studying A^A^ correlations. It is clear, from Figs. 121 and El that such a kinematics should 
be located around the two-nucleon correlation region, in order to exclude the influence of 
the final state interaction between the spectator nucleons. This requires high values of 
the missing energy and momentum of the active nucleon (see Fig. 0]). Another important 
condition is that the range of |q| and go should not be too far from the quasi-elastic peak, 
where x ~ 1. In this case the corrections from the off mass shell effects and meson production 
are minimized, and only the final state interaction of the hit nucleon with the spectators 
becomes relevant. However, it should always be kept in mind, that if a region exists where 
the interaction in the spectator pair (the A — 1 system in case of complex nuclei) can be 
neglected, this is no guarantee that the interaction of the struck nucleon with the nucleons 
of the spectator pair (the A — 1 system), can be neglected as well. It is clear therefore that 
one has to go beyond the PWIA, which is precisely the aim of the present paper. The effects 
of the full FSI on the process ^ife(e, e'2p)n have been recently investigated, treating the FSI 
within the GA 40j. In the present paper we will investigate the same topic in A{e,e'p)X 
process off and ^He within both the GA and the GEA. 



IV. THE FULL FINAL STATE INTERACTION WITHIN THE GENERALIZED 
EIKONAL APPROXIMATION 

Let us consider the interaction of the incoming virtual photon, 7*, with a bound nucleon 
(the active nucleon) of low virtuality (p^ ~ Mf^) at a kinematics not very different from 
the quasi-elastic one, i.e. corresponding to x ~ 1. In the quasi-elastic kinematics, the 
virtuality of the struck nucleon after 7*-absorption is also rather low and, provided pi is 
sufficiently high, nucleon rescattering with the "spectator" A — 1 can be described to a 
large extent in terms of multiple elastic scattering processes in the forward direction (in 
the system of reference where the target nucleon is at rest). These rescattering processes 
are diagrammatically depicted in Fig. El where, as in the rest of this paper, the internal 
and intermediate state momenta are denoted by fcj's and the final state momenta by pj's. 
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The diagrams essentially describe the process of multiple scattering in the most general case, 
within the assumption that all intermediate nucleons are on-shell. The low virtuality (before 
and after 7* absorption) of the active nucleon, coupled with the forward propagation, allows 
one to simplify the description of the final state interaction, which can be treated within 
the eikonal approximation. Before illustrating in detail the approach we have used to treat 
FSI in (e, e'p) reactions, we would like to discuss an important related issue (see also |4l|), 
namely the validity of the factorization approximation, frequently used in calculations at 
high Q^, and consisting in factorizing the {e,e'p) cross section into an e.m. and a nuclear 
parts, in spite of the fact that factorization, holding exactly in PWIA, is violated when FSI 
is taken into account. In the next Section the factorization approximation will be discussed 
within the GA and the GEA. 



A. The FSI in A{e, e'p)B processes within a diagrammatic approach 

Most of the problems one faces when trying to develop a fully covariant treatment of 
FSI, arise because of the hadrons' spins. Therefore, let us rewrite the hadronic tensor (Eq. 
0) in the following, fully equivalent form, which however exhibits explicitly the dependence 
upon the spin quantum numbers 

^^^ = j\rT. E Tl{MA.MA-l,S,)T,{MA.MA-l.s{){2T^Y5^'\PA + q-PA-l-pi), 

(13) 

where is a short-hand notation for the transition matrix element 

T^{Ma,Ma-i,si) = {aA-iPA-iEi_^,sipi\j;^{0)\aAPA). (14) 

The basic assumption underlying the eikonal diagrammatic treatment of FSI at high is 
that the transition matrix element for a nucleus A can be written in the following form 

A-l 

T^{Ma,Ma-i,si) = ^T^")(A^a,-Ma-i,si), (15) 

ra=0 

where the superscript (n) corresponds to the order of rescattering of the struck particle with 
the A — 1 nucleons (the "spectator" nucleons), namely Tjf^ corresponds to the PWIA (no 
rescattering), T^^^ to the single rescattering of the struck nucleons with the spectator ones, 
Tjf^ to double rescattering, and so on. Such an approach is expected to be vahd either at 

10 



high energies, when particles propagate mostly in the forward direction along the direction 
of the three-momentum transfer q, or when the momentum of the struck nucleon pi relative 
to A — 1 is sufficiently high; in both cases the eikonal approximation could be applied. 
The calculation of the rescattering part of in terms of Feynman diagrams appears in 
principle to be a prohibitive relativistic task due, as previously stressed, to the treatment 
of the spin. A relevant simplification occurs if the cross section factorizes into the e.m. and 
the nuclear parts and, as a matter of fact, many calculations performed within the eikonal 
approximation treatment of the FSI, simply assume factorization. Let us try to analyze the 
limits of validity of such an assumption, and to this end let us consider the deuteron. In this 
case the Feynman diagrams describing rescattering are given in Fig. El and the corresponding 
matrix element is 

T,{M2, si, S2) = T^'\M2, s,, S2) + T«(A<2, si, S2), (16) 

where Ai2, Si and S2 are the spin projection of the deuteron, and of nucleon "1" (the active 
nucleon) and nucleon "2" (the spectator nucleon) in the final state. Eq. (fTT)|l obviously 
states that in the deuteron the interaction between the struck and the spectator nucleon can 
occur only via single rescattering. The cross section of the process is given by 



dE'dn' '^^(27r)32Ei (27r)32E2^ 



where Ei = y pf + M^, and the hadronic tensor is as follows 

L?^ = ^\ E TliM2,s„S2)UM2,s„S2){27crS^'HPD + q-Pi-P2). (18) 

^ M2,Sl,S2 

Let us now obtain the factorization of the cross section, expressed in terms of the hadronic 
tensor (jl8|) . within a fully covariant approach. 



1. The PWIA and the factorization of the cross section. 



The PWIA for the process ^-ff(e, e'p)n within a covariant Feynman diagram approach has 
been considered by various authors (see e.g. [42, 43, 44 1). The matrix element = T^^-*, in 
such a case, has the following form 

T;,'\M2,SuS2) = 

^ E ^^(^''Pi' ^1' ^1) s~i)^D'ih, k2){k2 + M^)t;(k2, 82)] , (19) 



2M, 



N 
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where 



(20) 



TJ^*^{Q'^,kf) is the e.m. vertex, and $;^^(A;i, ^2) the covariant deuteron amphtude corre- 
sponding to the D —>■ (pn) vertex. The exphcit form of the amphtude ^^^{ki, ^2) depends 
upon the specific covariant model used to describe the deuteron and could be found else- 
where (see e.g., 42, Here, without loss of generality, we will use the Bethe-Salpeter 
(BS) formalism according to Refs. Q] and [4^ . 

When Eq. (fTUI) is placed in Eq. (fTH|) . the e.m. and nuclear parts gets coupled by 
the summation over the intermediate spins si and si'. However, it can be shown (see 
Appendix El that the square of the expression in brackets in Eq.()19p after summation over 
A^2 and S2 yields a 6 function 6~^-' , i.e. becomes diagonal in si, this leads to the decoupling 
between the e.m and the nuclear parts in Eq. p9|l . with the resulting hadronic tensor given 

by 

1 1 



-D 



2Md3 



J2 Tl{M2,Si,S2)T,{M2,SuS2) 



M2,Sl,S2 



2Mo (2E^,2E^,Lf^ {Q^p,, k,)) n^dkil) (27r)^5(^) (p^ + g _ _ p^) . (21) 



In Eq. fl^ . ud is the deuteron momentum distribution given by 

^i^(lkil) = ^ E \\u{^i,Si)<l>^'{ki,k2){k2 + M^)v{k2,S2) 



M2,S1,S2 



3 E \{h,S2\K'{k^: 



^(4(|ki|)+«l,(|ki|)), (22) 



M2,Sl,S2 

where the (covariant) deuteron wave function has been cast in a form similar to the non 
relativistic one with the scalar parts of the wave function, Ui(|k|)'s, related to the corre- 
sponding vertex functions, Giikf, k^ = M^), by a well known definition (see Eqs. (B8) and 
(B9)) leading to 



«L(|k|) 



2Mn- 



-GL {\Kho = MD-E^) 
k^-Ml 



(23) 



Placing Eq. ^21^ in Eq. p7|l the well known factorized form for the cross section is obtained. 



VIZ. 



a 



dE'dn'dpr. 



K{Q\ X, p„) a'^ {Q\ p„)nz5(|ki|))5(go + - E^,+^ - E^J. (24) 
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We reiterate that factorization has been obtained because the sum over S2 and A^2 in (fTH)) 
leads to the appearance of a delta function 5s', s^, which means, in turn, that the square of 
T<». becomes diagonal . TMs particular (exact) result . part ot a mote general assertiou 
that within the PWIA the nuclear Spectral Function is always diagonal in spins 35j . Let us 
now consider FSI; in this case the tensor (jl8|) is off-diagonal in spins and factorization does 
not occur. However, we will show that under certain kinematical conditions, satisfied to a 
large extent by the GA and GEA, factorization can be recovered. 



2. FSI: the single scattering contribution and factorization of the cross section. 

Let us compute the second diagram of Fig. E| To this end, we introduce a two-nucleon 
scattering operator T in terms of which the elastic scattering amplitude /^^, describing the 
elastic scattering of two on-shell nucleons, will be defined as follows 



/s-iS2;si,s2(Pi> P2; ki, ks) = m(pi, si)m(p2, S2) T M(ki, Si)u{k2, S2), 



(25) 



which is obviously the free A^A^ scattering amplitude; for a bound nucleon one has in prin- 
ciple to consider off-shell effects but in the GEA no virtuality is considered; this could be 
done for example by the approach of Ref. Q|, by introducing cut-off form factors in the 
corresponding nucleon lines, which formally leads to two Feynman diagrams with different 
"nucleonic" masses. In presence of FSI, the transition matrix element is 



5 S2) 



(26) 



with T^o) given agam by Eq. and T« given by the following form (note that henceforth 
we always have ki = — k2) 

d'k2 /.W,;.,.,(Pl,P2,k;,k2) 



T«(-^2,Sl,32) 



E 



-51^1 ^2 



X 



X 



u 



(k;, si')n*^(Q', A;f )M(ki, si)l [n(ki, si)$^^(A;i, k2)v{\i2. h 



(27) 



The full matrix element will therefore be 



T^{M2,Si,S2) 



2M, 



TV 



E W'' Pi' P^' ^1) s\)<^'^'{k^, k2){k2 + M^)t;(k2, S2] 



+ 
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+ 



X 



2M, 



N 



E 



SlSl S2 



i(27r) 



-M% + ie 



X 



(28) 



When Eq. (j2Hl) is placed into Eq. (fTHj) . the resulting equation is not diagonal in the 
spin quantum numbers and factorization does not hold. Let us however consider the basic 
assumptions underlying the eikonal multiple scattering approach, viz.: 

1. the momentum transfer k, in the elastic rescattering is small and mostly transverse i.e. 



K = Pl - k'l = k2 - P2 ^ k2_L - P2_L = 



(29) 



2. the spin-flip part of the A^A^ amplitude is very small, which means that, taking into 
account Point 1, one can write 



fsrsr,sus,iPl, P2, K, k2) ~ Ssy,s^Ss,,s2f 



(30) 



which is realized either at high values of the three-momentum transfer q, or at high 
values of the momentum pi of the struck nucleon relative to the A — 1 spectator 
nucleons. 

If the above conditions are satisfied, Eq. ()27|) assumes the following form (cf. Appendix B) 

1 r d% /^^(a^±) 



2MmJ i{2TTY - Ml + le 



(31) 



and one can write 



2M, 



N 



E'^'^(<5^Pm,Pl,S~l,Sl 



X 



+ 



M(ki, Si)$^^(A;i, k2){k2 + Mjv)T;(k2, S2 



+ 



u{k,,Sl)<^D'{h,k2)v{k2,S2] 



(32) 



i{2ny k[^ -M% + ie 

It can be seen that r^°) and T^^) in Eq. ^ have very similar structures, except that in 
T^^) the vector k2 is now an integration variable, since k2 7^ P2. When Eq. (fTHj) is evaluated, 
with given by Eq. and assuming soft A^A^ rescattering ( low values of n±), the main 
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contribution in the integral over k2 results from the region where k2 ~ P2 and this, in turn, 
originates again a delta function 6sj^si (^^e Appendix |B)) and the hadronic tensor becomes 

= E Tl{M2,s„S2)-n{M2,s„s,){27r)'6^'^{Po + q-p,-p2)^ 



M2,Sl,S2 

I E [jfKQ^ Pm, Pi, S2, S^) ■ jf{Q\ p„, pi, i-2, Si) 



2 - 



X 



X 



+ 



E 



u(k2, Si)<l>^^(A;i, A;2)(A;2 + M^)t;(k2, ^2) 



k2=P2 



+ 



M(k2,5i)$^^(A;i,fc2)T;(k2,S2)] 



i(27r)4 A;;2 - + ie 
X (27r)45(^)(PB + g-pi-p2) 



X 



(33) 



and the factorization of the e.m. and the nuclear parts is recovered. Eq. (j33|l could be 
written in a more familiar form if one integrates over k2o by taking into account the pole in 
the amplitude <^^'^{ki,k2) {k2o = -Ek) and neglecting the pole from the active propagator, 
which is located at large values of k2o and does not contribute to the integral. Using 
one obtains {J d^k/[i {27f)^] ^ J d^k/[{2E^{27if]) 



2Ep,2Ep„L;l(Q^p^,p 



X Y: (^l,S2K^(k2 

X (27r)V^)(Pz5 + g-pi-p2 



k2=P2 



+ 



X 



d^k-? 



2E^,{2T:f k[^ - Ml + le 



■Si, 521^/1) {iS~2 



X 



(34) 



By placing the above equation in Eq. (fT7|) . one obtains 
d^a 



K{Q\ X, p^ )a^^(g^ vM'^Vn.m^ + Md- ^p,+q - ^pj, (35) 



dE'dQ'dpn, 

where the Distorted Momentum Distribution n^^^ is 



E 



M2,S1,S2 



(.l,.2K^(k2)) 



k2=P2 



d^k. 



(Si,S2|vI/^^(pJ) 



2Ek2(27r)3A;i2-M2, + ^£ 



1 



J2 \'rf,'\M2,Si,S2)+r^'\M2,Si,S2) 



M2,S1,S2 

and the quantities 



(36) 



ri°)(-M2,Sl,S2) = (si,S2|^^^(k2 



(37) 



15 



and 

can be called the reduced (Lorentz index independent) amplitudes; in the above equations 
is the deuteron wave function and the spin wave function refers to the two particles in 
the continuum. 

To sum up we have shown that: 

1. the cross section which includes FSI factorizes provided: i) the spin flip part of the NN 
scattering amplitude can be disregarded, which is consistent with the high energies we 
are considering, and ii) the momentum transfer k, in the A^A^ rescattering is small and 
transverse, so that in the integral ()31|) one has k2 ~ P2 or, equivalently, k2 ~ Pm', this 
is a reasonable approximation, thanks to the behaviour of the elastic A^A^ scattering 
amplitude, which is sharply peaked in the forward direction; 

2. in the eikonal approximation and neglecting the spin dependence (spin-flip part) of 
the A^A^-amplitude, the FSI is not affected by the spin structure of the wave functions 
of the deuteron and the two-body final state. This means that in computing the 
Feynman diagrams, the intermediate spin algebra can be disregarded, and only the 
scalar part of the corresponding vertex functions can be considered, using Eq. (j23|l to 
define the scalar parts of the wave functions. Then the resulting amplitude has to be 
merely sandwiched between the spin functions of initial and final particles. 

These conclusions can be generalized to a nucleus A, for which the cross section of the 
process A{e, e'p){A — 1) is given by the following expression: 

^E,^^,^^^ = ^iQ'^ ^' P-) ^'""(Q'^ PjPP'iPm, (39) 
where P^^^Pm, E^) is the Distorted Spectral Function 



PP\Pm, Em) , 2 -H ^ ^ 

[Zn) ZJa + 1 j Ma,Ma-i,si 



2 

n=0 



X 



X 6i^Em-iE{_, + Emin)) (40) 

and n denotes the order of rescattering. In what follows the distorted momentum distribu- 
tions for the Deuteron (Eq. (IHSl)) and the distorted Spectral Function for ^He (Eq. (jlUj) ) 
will be calculated within the GA and GEA. 
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B. The process ^i?(e, e'p)n within the GA and GEA. 



Let us now calculate the reduced amplitude Tj^^ in the process ^i? (e, e'p)n, taking FSI 
into account by the GEA. This amounts to replace the energy denominator in Eq. (j38p by 
its generalized eikonal approximation. To this end, we will consider both the "canonical" 
case, when the value of the 3- momenta transfers |q| is so high that q pi, with the 2;-axis 
naturally directed along q, as well as the case of smaller values of q, but high values of pi, 
when q and pi may point to different directions, in which case the z-axis is oriented along 
Pi 

Remembering that k, = pi — k[ = k2 — P2, the energy denominator can be written as 
follows 



2|pi| 2|pii; 

2|Pi| - ^^^^ii^^«:ol ^ 2|pi| (k, + A,) , (41) 



2|pi 



where 



and the relation 



^. = ^^^^^^^{E^-\Ea\) (42) 
2|pi| 



k^ = E^,-E^,^^7^-{E^-\Ea\) (43) 



resulting from energy conservation go + = Ep^ + Ep^ has been used. 

By changing the normalization of the A^A^ amplitude from the covariant one to the non 
relativistic analogue {Ep ~ M^r), one has 



4Ep|pi| ^ 4M^|p 



L| = a^^(At^) = t J rf^be*'^^^ r(b) (44) 



and T^^^ becomes 



r^\M2,s„S2) = I yfl^a^Vx) J. ^. (^i,^2K-(p)). (45) 

Using 

^ ^2 = j ^(-2) e"^^''^^'^'"dz (46) 
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we obtain, in coordinate space, 

ri'\M2, S2) + tI,'\M2, S2) = { si, S2 (1 - ^(^)e^^^^r(b)) e-P-'-lM/^^W) . (47) 
As a result, the cross section will read as follows 

-- K{Q\x,Pm) a'''{Q\pm)n'o'''{Pm)S{MD + ly - Ep, - EpJ (48) 



dE'dQ'dpm 

with the distorted momentum distributions n^^^ defined by 

2 



nn' ip^m) = E / drxi,K'\r)Si'\r) exp(-^p^r) 

'^^^^1 M2,S23 '' 



(49) 



where S^^\t), which describes the final state interaction between the hit nucleon and the 
spectator, is 

Sl^\r) = 1 - e{zy^^'T{h) (50) 

with r = (b, 2;). In the above formulae the z-axis is along pi; it should be pointed out, 
however, that at large values of the momentum transfer, the hit nucleon propagates almost 
along q so that by choosing the ^-axis along the three-momentum transfer and neglecting 
the virtuality of the struck nucleon before and after interaction, one can write 

k'^ - Ml = ih + qf -Ml^ 2|q| (/«, + A,) , (51) 

where 

A, = (52) 

|q| 

It can be seen that the FSI factor (j50|l in the GEA differs from the one of the standard GA 
{4^, 42, 3|, simply by the additional factor e*^^^. It should be pointed out that whereas the 



well known factor 9{z) |42l,|4S| originates from the non relativistic reduction of the covariant 
Feynman diagrams and guarantees the correct time ordering of the rescattering processes, 
the quantity is of a pure nuclear structure origin and, as it can be seen from Eq. (|Tfj). 
represents a correction to the parallel component of the missing momentum. Therefore the 
corrections from A^ are expected to be important in parallel kinematics at |pz| ~ A^. As we 
shall see from the results of our calculations performed in perpendicular kinematics in the 
range |pm| < QOO MeV/c and Em < 100 MeV, one always has [A^l <^ |p|^ with go/|q| — 1, 
so that Aj, is always very small. We can therefore anticipate that effects of A^ on the 
experimental data we have considered is also very small. 
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C. The processes ^H{e,e'p)'^H and ^H{e,e'p){np) within the GA and GEA. 



Let us now consider the three-body system. The distorted Spectral Function is given by 
Eq. PHl 
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/ M3,M2, SI 



n=0 



X 



X 5 (Em - + Emin)) : (53) 



where the magnetic quantum number refers either to the deuteron or to the two nu- 
cleon in the continuum, depending upon the break-up channel we are considering {Emin = 
— E2 {Emin = E3) for the two-body (three-body) break-up channel). The diagrams rep- 
resenting the rescattering processes are shown in Fig. 6. The evaluation of these diagrams 
follows the standard procedure adopted for the deuteron. Let us illustrate it in the case 
of the 3bbu considering, for ease of presentation, the single scattering diagram of Fig. 6 
b). After integration over ^20 and k^Q in the corresponding poles of the propagators of the 
spectators {k2o = E^^ and ^30 = -Eks), we obtain 



2Ek,(27r)3 2Ek3(27r)3 

<^ge^l(23)(fel, ^2, fc3,gl,g2,g3) fNNjPl - k[) ^'(23)^/ (^2? ^3, -§2, S3) 

{kl-M%) kf-M% {k'i-Ml) 

where the overlaps of the vertex functions Gi are 



(54) 



GHe^l(2'i){kl, k2, fcs, Si, S2, S3) = (ki. Si, k2, S2, ks, S3\GHe^l{23){M3, P3)); (55) 

G(23)~.f{k2, k3, S2, S3) = (k2, S2, k3, S3 1 ^(23)^/ ( A^23, 5'23, P2, E[)); (56) 

The vertex functions Gi are replaced by the non relativistic overlap functions according 
to the general convention (we omit for ease of presentation the proper normalization factors) 

/c c c Ixij^Hh k k n ~ GHe^l(23){kl,k2,k3,Si,S2,S3) , . 

\Si,S2,S3|W^g (_ki,k2,K3j; (A?2~riVf^) ^ ^ 

and, using the completeness relation when summing over S2 and S3, one gets 
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X ^(,3)(k3, k- ; 5,3) j;^^^^ .^^ (^iK^(ki, k2, ka)). (58) 

Following the procedure adopted for the deuteron, we obtain 



where 



23 J 



d^K „ , , „ . /iviv('«)/4M;v|Pi 



(27r) 



3^(23) (k3,k^;-S23)^ 



(k^ + + ie) 



(5lK/(ki,k2,k3)), 



(59) 



-5'ki+q + -fi'pi 

2|Pi| 



(-Sm — £^3) ■ 



(60) 



Including also the 2buu channel, we can write, in coordinate space. 



P|f (p^, E^) = Pf/^(p^, E^) + Pr^(p^, E^), 



(61) 



where 



p//^(p^, ii;^) = <;r^^(p^)5(E„ - (^3 - ^2)) 



gr 



(62) 



with 



</'(Pm) 



1 1 



(27r)3 2 



E 1/ e^^P-xt,,^^^^(r)5i^^(p,rX/(p,r)dpdi 



(63) 



and 



n™(P™>-B".) 



X 



X 5 LB™, - 



— Er 



N 



(64) 



The FSI factor ^S^*^^ describes the single and double rescattering of nucleon "1" with the 
spectators " 2" and " 3" , and has the following form 
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S^'iP, r) = 5(^f (p, r) + (p, r) (65) 
with the single scattering contribution 5^^^^ given by 

(p, r) = 1 - ^ 9iz, - zi)e^^^(^-^^)r(bi - b,) (66) 

i=2 

and the double scattering contribution by (Ref. 13, [3|). 



+ eiz3 - Zi)e{z2 - ;,3)e-*^2(^3-.i)e-(A.-A.)(.,-.0] X r(bi - b2)r(bi - bg), (67) 

where Aj = (go/|q|)(-E'pi - E^,) and is given by Eq. 

i 

When A2 = 0, the familiar form for S^^^ is obtained, namely 

S^'\p, r) = n [1 - e{z, - z,) r(b, - bi)] , (68) 

i=2 

and when F = 0, the distorted Spectral Function (|UT|) transforms into the usual Spectral 
Function Q. 

Using Eq. (jHH), the cross section of the process ^He{e,e'p)X {X = D or (np)) assumes 
the following form 

^E'dn'dprr, = K{Q\x,pJa^''{Q\pm)P^riPm,Ej. (69) 

V. RESULTS OF THE CALCULATIONS 

We have used Eqs. (j^OI), and (jFTfjl to calculate the cross sections of the processes 

'^He{e, e'p)n, ^He{e, e'pYH and ^He{e, e'p){np). All calculations have been performed using 
the following well known parametrization of the profile function F(b) 

F(b) = <M}^^Mle-^y2t^,^ (70) 

where cr^^ is the total A^A^ cross section, a^N the ratio of the real to imaginary part of 
the forward A^A^ amplitude, and bo the slope of the differential elastic A^A^ cross section. 
The values of the energy dependent quantities a^j^% and a^N have been taken from Ref. 
3|. For the electron-nucleon cross section a'^^{Q'^, pm) we used the De Forest o'^^{Q'^, Pm) 
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cross section j35|. All two- and three-body wave functions are direct solutions of the non 
relativistic Schrodinger equation, therefore our calculations are fully parameter free. 

Calculations have been performed in PWIA and including the full rescattering within 
the GA and GEA by evaluating the Feynman diagrams shown in Figs. 5 and 6. It should 
be pointed out that, apart from minor differences (e.g. the structure of for complex 
nuclei) which do not affect the numerical results, our GEA is essentially the same as the one 



developed in 
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A. The process ^if (e, e'p)n 

Our results for the process ^-ff (e, e'p)n are compared in Figs. [7llHland|niwith three different 
sets of experimental data, covering different kinematical ranges, namely the experimental 

isl. 



data from NIKHEF 17], SLAG Hg, and Jlab 
e expe 

1 y 



The relevant kinematical variables in the 
three experiments are as follows: i) 0.1 < < 0.3, 0.3 < x < O.eJlJ; ii) 1.2 < < 6.8, 



X ~ 1 [13; iii) - 0.665 {GeV/cf , |q| - 0.7 GeV/c, x - 0.96 In Figs. Elandinithe 
theoretical cross section corresponding to Eq. (jlH)) . namely 

= free /C(Q^ X, p^) at^,{Q', p„)ng^^(p„) (71) 



dE'dQ'dn^ 



is compared with the corresponding data, whereas in Fig. |Hlwe compare, as in Ref . |l8l| . the 
effective momentum distributions Neff{pm) (or reduced cross section) defined by |l8| 

Pm 

where in Eqs. (\71\i and (f?^ free and /C are kinematical factors which arise from the inte- 
gration over dTpj^. 

The results presented in Figs. 7, 8 and 9 exhibit a general satisfactory agreement between 
theoretical calculations and experimental data, particularly in view of the wide range of 
kinematics covered by the data we have considered. Figs. 8 and 9 show however that 
quantitative disagreements with data exist in some regions. Particularly worth being noted 
is the disagreement in the region around |pm| — 0.25 GeV/c appearing in Fig. 8. We did not 
try to remove such a disagreement by adjusting the quantities entering the profile function 
(I48|l . but it turns out that n^^^ in the region around |pm| — 0.25 GeV/c is rather sensitive 
to the value of a^N- The NIKHEF kinematics deserves a particular comment. As a matter 
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of fact, the four momentum transfer in this experiment is rather low and one might rise 
doubts as to the vahdity of the eikonal approximation. In this respect, it should however 
be pointed out, that what really matters in GA and GEA is the relative three- momentum 
of the hit nucleon with respect to the A — 1 system; in the NIKHEF experiment, due to 
the large value of the energy transfer, the three momentum transfer is also large, and 7* 
absorption occurs on a proton moving along q, with the recoiling neutron moving with low 
momentum against q; the resulting proton-neutron relative momentum is of the order of few 
hundreds MeV/c, which, though representing the lower limit for the validity of the eikonal 
approximation, still appears, according to our results, to be suitable for the application of 
the GEA: as a matter of fact our results appear to be in reasonable agreement with the ones 
obtained within ap pro aches which are better justified at low energies, like, e.g., the ones 
presented, in Refs. H H 52 , 3, 3 • 

Let us conclude this Section by stressing that, as far as the effects of MEG and A isobar 
excitations are concerned, these have been found to be very small (~ 5 — 6%) in the SLAG 

n n 

kinematics ( see with the results from exhibiting the same trend also in the Jlab 

kinematics. Eventually, we should remark that the results of the GA and GEA differ by 
only few percent and cannot be distinguished in the Figures. 



B. The processes ^ He{e,e'p)'^H and ^He{e,e'p){np). 

Galculations for the three-body systems are very involved, mainly because of the com- 
plex structure of the wave function of Ref. which is given in a mixed [Lp, X, j23, S23) 
representation, including angular momentum values up to Lp = 7 and j23 = 8 (a total of 
58 configurations with different combinations of {Lp, X, j23, S23) quantum numbers). Gorre- 
spondingly, the wave function of the spectators (the deuteron or the continuum two-nucleon 
states) is given in a JLS'-scheme (see Appendix We would like to stress, that no approx- 
imations have been made in the evaluation of the single and double scattering contributions 
to the FSI: proper intrinsic coordinates have been used and the energy dependence of the 
profile function has been taken into account in the properly chosen GM system of the in- 
teracting pair. The Feynman diagrams which have to be evaluated, both for the 2bbu and 
3bbu channels are shown in Fig. 6. 
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1. The two-body break-up channel ^ He{e, e' pf' H . 



The 2bbu channel cross section 



1- (|nni), is compared in Fig. 10 with 
21l. The relevant kinematical variables 



obtained from Eq. with n^/'^(pm) given by 

recent experimental data from Jlab Collaboration 
in the experiment are |q| = 1.5 GeV/c, go = 0.84 GeV , = 1.55 {GeV/cY, and x ~ 1. 
The cross section is presented as a function of the missing momentum \pm\ (which, for the 
^i7e(e, e'p)Z^-process, exactly coincides with the final deuteron momentum). In PWIA the 
cross section is directly proportional to Ugr (Eq. ()lip) shown in the left panel of Fig. 01 
It can be seen that up to |pm| ~ 400 MeV/c, the PWIA and FSI results are almost the 
same and fairly well agree with the experimental data, which means, in turn, that the 
2bbu ^He{e,e'p)'^H does provide information on Ugr] on the contrary, at larger values of 
|Pm| > 400 MeV/c the PWIA appreciably underestimates the experimental data. It is very 
gratifying to see that when FSI is taken into account, the disagreement is fully removed and 
an overall very good agreement between theoretical predictions and experimental data is 
obtained. It should be pointed out that the experimental data shown in Fig. correspond 
to the perpendicular kinematics, when the deuteron momentum (the missing momentum) is 
always almost perpendicular to the momentum transfer q; in such a kinematics the effects 
from FSI are maximiz ed, whereas in the so called parallel kinematics, they are minimized 
(see, e.g. j^, 5g). The kinematics therefore reflects itself in the relevance of the 

calculated FSI; as a matter of fact, we have found that the effects of the FSI calculated either 
within the GA or GEA approximations, differ only by a few percent, which was expected in 
view of the observation that the factor (Eqs. (jUUI) or affects only the longitudinal 

component of Pm and therefore has minor effects on the data we have considered. The effects 
of MEC and A isobar contributions have been estimated in and found negligible up to 
about Pm ^ QOOMeV/c. 
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2. The three-body break-up channel ^He{e,e'p){np). 



From Eq. (j6yj), we obtain the cross section for the 3bbu in the following form 

^^'^ '--^ ^ eN//^2 „ \r>FSI 



KshbuiQ, X, Pm) P™)Pe/ (Pm, Ej, (74) 



dE'dn'dnp.dE^ 

where Pf^^^ {pm, Em) is given by Eq. (jMj) . We have calculated Eq. (fTIj) in correspondence 
of two different kinematical ranges: the one from Ref. [2^ and the one corresponding to 
the experimental data from Jlab Contrary to the 2bbu channel, the 3bbu cross section 
depends upon an extra kinematical variable, the removal energy Em, and corresponds to 
the process in which three particles interact in the continuum. We have considered three 
different theoretical approaches, namely: 

1. the Plane Wave Approximation (PWA), when FSI effects are completely ignored , i. 
e. the three particles in the continuum are described by plane waves; 

2. the Plane Wave Impulse Approximation (PWIA), in which the struck nucleon is de- 
scribed in the continuum by a plane wave and the spectator pair is described by the 
continuum solution of the Schrodinger equation (obviously, in the case of the deuteron 
the PWIA coincides with the PWA); 

3. the full FSI, when the struck nucleon interacts in the continuum with the nucleons of 
the spectator pair via the standard GA or the more refined GEA. 

In Fig. p the results of our calculations are compared with the experimental data 
from Ref. j20|. In the experiment, which corresponds to a relatively low beam energy 
{E = 0.560 GeV), the scattering angle {9e = 25°) and the energy transfer (go = 0.32GeV^) 
were kept constant, and protons with different values of the missing momentum and en- 
ergy were detected in correspondence of several values of the proton emission angle 6'pj, 
viz O-p^ = 45", 60°, 90.5°, 112° and 142.5°. The kinematics is far from the quasi elas- 
tic peak (x ~ 0.1) and the values of the four- and three- momentum transfers are low 
(Q2 ~ 0.03 {GeV/c)'^ and |q| = 0.28 GeV/c). At first glance this would invalidate the 
use of the eikonal approximation; however, a detailed analysis of the kinematics, shows that 
the value of both pi and Pm are rather large (400 — 600 MeV/c), and so is the value of the 
angle between them (6r~~ ~ 150°); thus the momentum of the struck nucleon relative to 
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the spectator pair is high enough to make the use of the eikonal approximation justified. 
Moreover, the values of the experimentally measured missing momenta and missing energy 
at each value of 6p^ , always cover the kinematical range where the condition for two nucleon 
correlation mechanism E^r, ~ — — holds; as a matter of fact, as it can be seen from Fig. 
11, the position of the bumps in the cross section are reasonably predicted by the PWA and 
PWIA. 

The results presented in Fig. El clearly show that with increasing missing momentum, 
the experimental peak moves to higher values of missing energy, in qualitative agreement 
with the two- nucleon correlation mechanism. More important, it can be seen that at the 
highest value of \pm\ {Op-^ = 112°) the effects of FSI, both in the spectator pair and between 
the struck nucleon and the spectator pair, is very small. The reason for such a behaviour 
is as follows: the kinematics of the experiment is not purely perpendicular: the relation 
between |pm_L| and |pm| is such that |pm±| ~ 2lP'"l' ^^^^ dominant role played by 
FSI in the purely perpendicular kinematics is decreased with increased values of 9p^ . 

In Fig. El our results are compared with the recent data from Jlab [22], where the 

cross section was measured at fixed values \pm\ vs. the missing energy Em- As in the 

case of the Saclay data previously analyzed, even in this case the cross section exhibits 

bumps approximately located at values of Em and |pm| satisfying the two- nucleon correlation 

p2 

mechanism relation {Em ~ — and in agreement with the behaviour of the Spectral 
Function (see Figs. ElandE]). However, contrary to the Saclay case, the PWIA dramatically 
underestimates the experimental data. This is clear evidence that the FSI between the struck 
nucleon and the nucleons of the spectator pair (Feynman diagrams b) and c) in Fig. ^ does 
play a relevant role, as the results of our calculations (the full fine in Fig. IT^ do indeed 
really show. Since, as already stressed, the Jlab experiment correspond to a perpendicular 
kinematics, this explains the larger effects of the FSI with respect to the Saclay experiment. 
The effects of the FSI calculated either within the GA or GEA approximations, differ only 
by a few percent, which was expected in view of the observation that the factor (Eqs. 
()fi()|) or ()52|) ) affects only the longitudinal component of Pm and therefore has minor effects 
on the data we have considered. 

There exist at present only two approaches to the calculations of the 2bbu and 3bbu 
channels at the Jlab kinematics: the one presented in this paper and the one by Laget 



reported in Refs. 



n 



22| and j23[. A comparison of the results of the two approaches 
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exhibits an encouraging agreement both in the 2bbu and 3bbu channels, with some minor 
differences which should most likely be ascribed to the different wave functions used in the 
two calculations. It is therefore gratifying to observe that different approximations to the 
treatment of FSl lead to very similar results. 

The effects of MEC and A, as previously pointed out, have not yet been considered in 
our approach; the calculation of Ref. Q], shows they reduce the cross section in the peak 
by about 10%, leaving the missing energy dependence and, consequently, our conclusions, 
practically unchanged. 



VI. SUMMARY AND CONCLUSIONS 

We have calculated the cross section of the processes ^H{e, e'p)n, ^iJe(e, e'p)D, and 
^He{e,e'p){np), using realistic wave functions for the ground state, which exhibits the very 
rich correlation structure generated by modern NN interactions; the FSI of the struck nucleon 



with t 



le spectators has been treated within the standard G 



auber eikonal approximation 



(GA) 3|, as well as with its generalized version (GEA) 0, 1^ The two approaches 
differ by a factor (Eqs. (jl^ and (IH2I)) which modifies (see Eq. (jHSJ) the FSI factor 
appearing in the standard GA (Eq. (jHHl))- This factor takes into account in the NN scat- 
tering amplitude the removal energy of the struck nucleon, or, equivalently, the excitation 
energy of the system A — 1. By properly choosing the 2;-axis (along q or pi), we were able 
to calculate FSI effects either in the case of large values of the three-momentum transfer q, 
or large values of the momentum of the struck nucleon pi relative to the A — 1 system; by 
this way calculations could be extended successfully even at relatively low values of Q^. As 
far as the three-body break-up channel in ^He is concerned, the FSI in the spectator pair 
was always calculated by the solution of the Schrodinger equation, whereas the interaction 
of the active, fast nucleon with the two nucleons of the spectator pair has been taken care 
of by the GA or GEA approximations. The method we have used is a very transparent 
one and fully parameter free: it is based upon Eqs. (jl^ . (jHH), and which only require 
the knowledge of the nuclear wave functions, since the FSI factor is fixed directly by NN 
scattering data. Of course with increasing A, the order of rescattering increases up to the 
[A — l)-th order; we have performed calculations in the three-body case exactly, and did not 
investigate the problem of the convergence of the multiple scattering series. This problem 
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is under investigation in the case of ^ife. Most of our calculations have been performed in 
kinematical conditions where the effects of MEC, A isobar creation, etc. are minimized, as 



confirmed by calculations performed, e.g., in Refs. j23|, |50|, |5l|, |52, |53|, |55|, |57|, |5^ As for the 

main results we have obtained, the following remarks are in order: 

1. the agreement between the results of our calculations and the experimental data for 
both the deuteron and '^iJe, is a very satisfactory one, particularly in view of the lack 
of any adjustable parameter in our approach; 

2. the effects of the FSI are such that they systematically bring theoretical calculations 
in better agreement with the experimental data. For some quantities, FSI simply 
improve the agreement between theory and experiment (cf. e.g. Figs. 7, 9 and 11), 
whereas for some other quantities, they play a dominant role (see e.g. Fig. 10 and 

12); 

3. a comparison of the PWA and the PWIA with the full FSI calculation, does show 
that proper kinematics conditions could be found corresponding to an overall very 
small effect of FSI, leaving thus room for the investigation of the details of the nuclear 
wave function; as a matter of fact, we always found that in the 3bbu channel in ^He, 
^He{e, e'p){np), the experimental values of Pm and corresponding to the maximum 
values of the cross section, satisfy to a large extent the relation predicted by the two- 
nucleon correlation mechanism j^Sl, namely Em — p^/4Mjv + E^ (cf. Fig. 3, right 
panel), with the full FSI mainly affecting only the magnitude of the cross section; 
thus, quasi elastic one-nucleon emission A{e,e'p)B processes at x ^ 1, together with 
processes at a; ~ 2, when the virtual photon is absorbed by a correlated two-nucleon 
" system" , would represent valuable tool for the investigation of correlations in nuclei; 

4. calculations of the 2bbu channel disintegration of '^He, i.e. the process ^ife(e, a'p^H, 
have already been performed Q| using realistic wave functions and taking exactly 
into account nucleon rescattering up to 3rd order, i.e by using the generalization of 
Eq. to the four-particle case, viz 

S^^^ = S^^^^K, ri2, r34) + 5(^^^(R, ri2, r34) + 5(^^^(R, r^, r^^), (75) 

where R, ri2, and are four-body Jacobi coordinates. Calculations for the 3bbu 
and 4bbu channels are in progress and will be reported elsewhere js^; they should 
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in principle yield results appreciably differing from the predictions based upon shell- 
model type four-body wave functions; 



5. our results for ^He generally agree with the ones obtained in Ref. j23[, so that it would 
appear that the problem of the treatment of FSI at high values of (or high pi) is 
under control; nevertheless, a systematic comparison of the various approaches would 
be highly desirable; 

6. we have given the criteria according to which at high energies the exclusive A{e, e'p)B 
cross section should factorize, and the similarity of our results withe the ones based 
upon a non factorized cross section [2^, confirm the validity of these criteria. 

7. eventually, it appears that in the kinematical range we have considered only minor 
numerical differences were found between the conventional Glauber-eikonal approach 
and its generalized extension; this does not mean at all that the same will hold in 
„tKe.Uine.a«ca>co„d.,o,.(seee.,Sya). 
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APPENDIX A: THE NUCLEAR WAVE FUNCTIONS 

In our calculations we have used two- and three-body wave functions corresponding to 
the AV18 potential 
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1. The ground state wave function of He 



For the ^He wave function we have adopted the correlated variational wave function 
by the Pisa group Q which is written in a mixed (Lp, X, J235 'S'23)-representation, where 
j23 and 5*23 are the total angular momentum and the total spin of the pair "23", X is an 
intermediate angular momentum resulting from the coupling j23 + Si and Lp is the radial 
angular momentum of the motion of the nucleon "1" relative to the pair "23". The explicit 
form of the wave function is 

^^:ip,r) = E E {XMx LpTUp |^^3) 023^23 \XMx YlMp) 

{a} {m} 

(^23/^23 S23Jy23\j23m23)Yl,,^,,{r)xS23U23R{a}{r, p)llf , (Al) 

2 2 

where {a} labels all possible configurations in ^He with quantum numbers Lp,X ^23,823, 
and T23 and (ZimiZ2?Ti2| luf^u) is a Clebsch-Gordan coefficient. The total isospin function 
is Jff = J2 { T23T23 |ti \\\ )Jt23T23V\tii where Jraaras and r^i^^ are the isospin functions 
of the pair and the nucleon, respectively. Obviously, because of Pauli principle and parity 
constraints, the allowed configurations in eq. ()A1|) are those that satisfy the following 
conditions: 

Lp + ^23 is even and I23 + S23 + T23 is odd (A2) 

The corresponding radial part of the wave function, _R{q,} (r, p) , has been obtained Q by 
a variational method using the AyiS potential including into the calculations values of 
Lp, Z12 = . . . 9 (a total of 58 different configurations Lp, X, J235 ^23, 'S'23 have been considered). 



2. The two-body continuum wave function \I'23(^^)' 

With the representation ()A1|) of the ^He wave function, it was convenient to adopt for the 
two-nucleon scattering state '^\.^{y) the spin-channel representation ^^21,^23^^) ■, characterized 
by the total (conserved in the scattering process) spin S'23 and its projection z/23. For spin 
5*23 = 1 one has 

vl>*^^,^(r) = 47r 5: 5: (/o/xo SfVf\JfMj )Y,„p„(t)i?J|,^,^(r)z'/3^iif'(i-):7T,3.23, (A3) 

where loilf = J/ ± 1, Jf. Note that the presence of tensor forces in the NN-potential 
leads to an admixture of partial waves with I = Jj — 1 and / = Jf + 1. This hinders the 
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use of real phase shifts for the asymptotic behaviour of the radial functions and, 
consequently, the Schrodinger equation cannot be solved in terms of real solutions. However, 
a unitary transformation V allows one to define new radial functions R = VR which are 
eigenfunctions of the scattering problem, i.e., solutions of the Schrodinger equation with the 
proper asymptotic behaviour. 



3. Wave function overlaps and the spectral function P{\hi\,E) oi^He. 



The Spectral function for the three-body break-up channel can be expressed in terms 
of the overlap between the three-body and two-body radial functions by substituting (jAlll - 
()A3|) into Eq. (fT^ . Using the orthogonality of the spherical harmonics Yim(t) and the 
completeness of the Clebsch-Gordan coefficients, one obtains that only diagonal ({«} = 
{aAf}) matrix elements contribute to the spectral function, viz. 



^e.(|kl|,E) = 2E E 

M:i af,Sf,Uf 



dH 



-ipki 



M, 



N 



{a} 



where fiso = 3(1) for the pair in the isosinglet (isotriplet) final state, jipipp) is the spherical 
Bessel functions, and the dimensionless overlap integrals O^^^^ (p) are defined as follows 

The normalization of the proton spectral function ()A4|) - ()A5|1 is 



d^kidEP{\ki\,E) 



0.15 for T23 = 
0.50 for T23 = 1 



(A6) 



so that the two-body break-up channel is normalized to ^1.35. Since the FSI factors 
S^^^ and S^^^ (Eqs. and (|UHj)) are no spherically symmetric, the distorted Spectral 

Function Pf^^^ {pm, E^) (EqEH) is not longer diagonal with respect to the {Lp, X, j23, S23) 
configurations. Except for parity constraints ()A2|) . any values of angular momenta of the 
pair in the final state contribute to Pf^^^{pm, E^)- 
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APPENDIX B: FACTORIZATION OF THE COVARIANT CROSS SECTION 



In this Appendix we will show, within a fully covariant approach, that under certain 
kinematical conditions the cross section for the process A{e, e'p)X process factorizes even in 
presence of FSI. We shall consider, to this end, the deuteron (D) treated within the Bethe- 
Salpeter (BS) formalism. As mentioned, the factorization depends upon the spin structure 
of the square of the matrix element n(ki, si)^^^{ki, k2){k2 + Mtv)^ (k2, S2) appearing in 
Eq. (fT^ or, in case of FSI, upon the structure of ^(ki, s'i)<l';^^(A;i, A;2)f (k2, S2) (cf Eq. 
(|2H1))- The relevant spin parts can be evaluated directly by using the explicit form of the 
Dirac spinors, u and f, and the explicit expressions for the amplitudes $;^^(A;i, ^2) (cf. 



Refs. j42. 
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1. The PWIA 

In Ref. the Feynman diagrams for the process D{e,e'p)X have been evaluated in- 
cluding all BS components. Here we re-calculate the diagrams for the ^S^^ and ^Df~^ 
components in a slightly different manner which will be useful when FSI effects are consid- 
ered. 

In PWIA the cross section reads as follows 

dE'dn' '^^(27r)32Ei(27r)32E2' ^ ' 

where l^'^ and are the leptonic and hadronic tensors, respectively, the latter being 

L'',. = 7^\ E 7;(-M2,si,S2)T,(>[2,si,S2)(2vr)V^)(Pz. + g-pi-p2)= (B2) 

= 7^1 E {M2\j;!\P2, S2,pu Si ) (pi, si,p2, S2I J,^|A<2 )(27r)35(^) {P^ + q - p, - P2) , 

^ M2,S1,S2 

where is the nucleon electromagnetic current operator. The amplitude could be 
written in the following form 

T,iM2, si, S2) = m(pi, si)r;^*^(g2, kl)<!>^'{ki, k2)s-\k2)v{p2, S2), (B3) 

where is a short-hand notation for the main BS amplitudes $3q++ and $3n++, corre- 
spending to L = and L = 2, respectively (see Refs. |42, |45[), ^2 = P2, S '^{k2) = k2 + m, 
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and TJ^* ^ {Q"^ , kf) is the electromagnetic eA^ vertex which, for an off-mass-shell nucleon, 
depends not only upon Q^, but upon kf ^ as well. 

By introducing between r^{Q'^,pi,ki) and $;^°(A;i, ^2) the complete set of the Dirac 
spinors 

1 



2M, 



N 



J2 [w(ki, si)u{ki, si) - v{ki, si)v(ki, si)] 



(B4) 



and bearing in mind that for the ^Sf^ and ^Df^ partial waves the second term in ()B4|) does 
not contribute, we obtain 
1 



T^iM2, si, S2] 



2Ma 



Y.jf{Q\pi,k,,~Si,Si) u{k,rsi)<^Diki,k2)S-\h)v{p2,S2) ,{B5) 



where Jf{Q\puh) = {p,, s,\rj:^{Q\kl)\k,,~Si). 
Let us evaluate Eq. ()B5|) for the D wave. One has 



n(ki, h)<^t'nMh, k2)S-\k2Hp2, S2) 



■^{kl - M^Mh, |k|) 2m {xi \[-{<Ti^) + 3K^)(na)} 

ki 



X. 



S2 / ) 



(B6) 



where n is a unit vector along k, i.e., n = — 

|k| 

expression for the cross section, one obtains 



|ki| 



. When Eq. ()B6|) is inserted in the 



\ E (xs-,|{-(^^-^^) + 3K-^^)(n^)} 



M2,S2 



XS2 [X. 



26. 



SlS' 



(B7) 



The last relation ensures factorization of the cross section; as a matter of fact, by performing 
the same procedure for the 5- wave, it easy to show that, thanks to Eq. ()B7j) . the cross section 
(Eq. (jBlll ) factorizes, assuming the form with no given by Eq. (j22I). In obtaining the 
above equations we expressed the BS amplitudes 4>L{ko, |k|) in terms of the BS vertices 
G3x^++{ko, |k|) and the radial functions ul, by the relations 



NNf 



{kl-M^^)Mko,\k\) 



NN^2EkG,jj++{ko, |k|) 
71 M^-2Ek ' 



(B8) 



where /cq = — — -Ek, and 



US{D) 



'^3g++(3g++)(fco, |k|)/(47r) 
v/2Md(Mb - 2Ek) 



(B9) 
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Note that in Eq. ()B9|) the normahzation of the wave function is chosen so as to correspond 
to the non relativistic deuteron wave function 

lj\k\'d\k\{uU\k\) + uU\k\))^l. (BIO) 

We reiterate that factorization in PWIA occurs because the sum over S2 and A^2 of the 
square of the matrix element in Eq. ()B7I ) becomes diagonal with respect to §1. 

When FSI is taken into account, instead of Eq. ()B7|1 . one obtains for the D-wave (for 
the 5*— wave the spin structure is trivial) 



I ^ [n(ki,si)$^^(fci,A;2)t;(k2,S2)]^[M(k;,si)<l>^Hfc;,A;^)t;(k^,S2) 
I E {Xs. + 3(n^^^)(no-)}| X..) {xs, |{-(^r^^) + 3(n'r^^)(nV)}| x.,) = 
I E {Xs. + 3(n|^^)(no-)} {-(^r-^^) + 3(n'^+^^)(nV)}| X.-,) (Bll) 



2 M2 

where n (n') is a unit vector along ki (k2). 

By taking into account the completeness of the polarization vectors the only spin 
dependence remaining in Eq. ()B11|) is contained in the term 

\ $:(n^'')(no-)(n'r^O(nV) = (nn') ( (nn') -^a[nx n']) , (B12) 

SO that in case of rescattering with low momentum transfer, when in the integral over k2 the 
main contribution comes from k2 ~ P2 , k2 ~ P2 , one has cr [n x n'] = and factorization 
is approximately recovered, with the S and D waves adding incoherently. 
Thus, to sum up, factorization is compatible with FSI if: 

1. The spin-flip part of the A^A^ amplitude should be very small, as it occurs when either 
the three-momentum transfer q, or the momentum |pi| are large; 

2. the momentum transfer k in the A^A^ rescattering has to be small small so that in 
the integral k2 ~ P2. This appears to be the case since the A^A^ amplitude is sharply 
peaked forward. 

3. The contribution from A^A" pair currents can be neglected, which is to a large extent 
legitimate due to the smallness of the P wave in the deuteron. 
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c) 

FIG. 2: Feynman diagrams for the process A{e, e'p){A— 1). (a) describes the Plane Wave Impulse 
Approximation (PWIA); (b) the single rescattering; (c) the full A — 1 rescattering. The four- 
momenta of particle i before and after rescattering are denoted by ki, k[, k'- , etc., respectively. 
The black oval spots denote the elastic nucleon-nucleon scattering matrix.. 
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FIG. 3: The proton Spectral Function of "^He (Eq.®). Left panel: Ugr (Eq. Hill) ) vs p = |ki|. 

Right panel: P^x (Eq. 1)12(1 ') vs the excitation energy of the two-nucleon system in the continuum 
e f 

Erei = -T-r- = E2 = E — E„iin, for various values of p = |ki|. The dot-dashed curves represent 
the Plane Wave Approximation (PWA), when the three particles in the continuum are described 
by plane waves, whereas the full curves correspond to the PWIA, when the interaction in the 
spectator proton-neutron pair is taken into account. The arrows indicate the position of the peak 
(~ p^ /AMn) predicted by the two-nucleon correlation model for the Spectral Function [s^ (three- 
body wave function from Ref. [2I, ^1/18 interaction 1^). 
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FIG. 4: The neutron Spectral Function of ^He (Eq. 1)12(1 ) vs p = |ki| and the excitation energy 
of the two-nucleon system in the continuum Erei = -r-r- = E2 = E — Emin- The dotted curves 
represent the PWA, when the three particles in the continuum are described by plane waves, 
whereas the full curves correspond to the PWIA, when the interaction in the spectator proton- 
proton pair is taken into account (three-body wave function from Ref. ^2], AVIS interaction ^). 
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FIG. 5: The Feynman diagrams for the process '^H{e,e'p)n representing the PWIA) (a)), and the 
single (6)) rescattering in the final state, f^^ denotes the elastic NN scattering amplitude. 



42 
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FIG. 6: The Feynman diagrams representing the PWIA (a)), the single (b)), and double (c)) 
rescattering in the processes ^He{e,e'p)D and ^He{e,e'p){np). In the former case the final two- 
nucleon state is a deuteron with momentum Pd = P2, whereas in the latter case the final state 
represents two free nucleons with momenta p2 and ps, with P2 = P2 + Ps- The trivial single and 
double rescattering diagrams with nucleons "2" and "3" interchanged are not drawn. The black 
oval spots denote the elastic nucleon-nucleon (NN) scattering matrix T (See Eq. (25)). 
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FIG. 7: The process ^H{e, e'p)n: the NIKHEF experimental data vs the missing momentum 
Pm = \Pm\ are compared with our theoretical calculations; the dotted line represents the PWIA, 
whereas the full line include the final state rescattering. The curves labelled LOl, L02, and L03, 
correspond to = 0.1, 0.2, and 0.3 (GeV/c)'^, respectively, and x ~ 0.3 — 0.6 (in this Figure and 
in Figs. 8-12, p' = |pi|). 
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FIG. 8: The process ^H{e, e'p)n: the Jlab experimental data 1^ {^eff defined by Eq. (f72|) ) vs 
the missing momentum pjji = | Pm 1 1 compared with our theoretical calculations . The dotted line 
represents the PWIA whereas the full line includes the final state rescattering. The experimental 
data correspond to the perpendicular kinematics, with ~ 0.665 (GeV/c)'^ , |q| ~ 0.7 GeV/c, 
and X ~ 0.96. 
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FIG. 9: The process H{e, e'p)n: the SLAG experimental data at x ~ 1 vs the missins; 
momentum pm = \Pm,\, compared with our theoretical calculations. The dotted lines represent the 
PWIA and the full lines include the final state rescattering. The positive and negative values of 
Pm correspond to values of the azimuthal angle (p = ir and 0, respectively. 
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!^IG. 10: The process ^He{e, e'p)D: the experimental data from JLab (JLab Experiment E-89-044 
2l|) vs pm = \pm\ compared, at = 1.55 (GeV/c)'^ and x = 1, with our theoretical results. The 

dashed line corresponds to the PWIA and the full line includes the full FSI calculated using Eq. 
the predictions by Eq.^^ (GEA) and Eq. ((HHl) (GA)^iffer by at most 4% and cannot be 

distinguished in the Figure (three-body wave function from [3], ^T^18 interaction jsl). 
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FIG. 11: The process^ He{e,e'p){np) : the experimental data from Saclay 
values of the proton emission angle 9p-^ (p' = |pi|), are compared with our theoretical results. The 
dotted lines correspond to the PWA, when the three nucleons in the final state are described by 
plane waves, the dashed lines correspond to the PWIA, when the interaction in the spectator 
neutron-proton pair is taken into account, and the full lines include the full FSI calculated using 
Eq. JMI); the predictions by Eq.^ (GEA) and Eq. ^ (GA), differ by at most 4% and cannot be 
distinguished in the Figure. Note that the values of the experimental pm and Em corresponding to 
the maxima of the cross section, satisfy to a large extent the relation predicted by the two-nucleon 
correlation mechanism js^, namely Em — Pm/^Mj^ (cf. Fig. 3, right panel), with the full FSI 
mainly affecting only the magnitude of the cross section ( three-body wave function from Q], ^FIS 
interaction jj]). 
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FIG. 12: The process ^ He{e,e' p){np)\ the same as m Fig. but for the Jlab experimental 

data j22i]) at = 1.55 {GeV/c)^ and x = 1. Note, that unhke what shown in Fig. 1111 here 

the differential cross section is plotted vs the excitation energy of the two-nucleon system in the 

f 

continuum i.e. E're/ = = £"2 = Em — E^. 

Mm 
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